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We use microscopic linear response theory to derive a set of equations that provide a complete 
description of coupled spin and charge diffusive transport in a two-dimensional electron gas (2DEG) 
with the Rashba spin-orbit (SO) interaction. These equations capture a number of interrelated 
effects including spin accumulation and diffusion, Dyakonov-Perel spin relaxation, magnetoelectric, 
and spin-galvanic effects. They can be used under very general circumstances to model transport 
experiments in 2DEG systems that involve either electrical or optical spin injection. We comment on 
the relationship between these equations and the exact spin and charge density operator equations 
of motion. As an example of the application of our equations, we consider a simple electrical 
spin injection experiment and show that a voltage will develop between two ferromagnetic contacts 
if a spin-polarized current is injected into a 2DEG, that depends on the relative magnetization 
orientation of the contacts. This voltage is present even when the separation between the contacts 
is larger than the spin diffusion length. 



I. INTRODUCTION 

Spintronics is an active field which studies processes 
that manipulate and probe the electronic spin degree- 
of-freedom, with the goal of identifying effects that can 
augment the orbital control and measurement procedures 
used in traditional electronics A Spin-related transport ef- 
fects in ferromagnetic metals are already used in current 
technology to provide the robust and responsive magnetic 
field sensors required by magnetic information storage 
systems. Hopes that spin-related transport effects might 
play a greater role in future information processing and 
storage technologies have motivated a growing body of 
research on the creation of spin-polarized carrier distribu- 
tion in semiconductors, either optically^ or by injection 
from other magnetic systems^^ Semiconductor quan- 
tum well electron gas systems are especially promising 
for spintronics because their intrinsic spin-orbit (SO) in- 
teractions are weak, implying long spin memory times, 
and because the Rashba SO interaction^ which enables 
electrical control of spin, can be tuned over a wide range 
by applying growth-direction gate potentials^ 

Although the study of spin dynamics in semiconduc- 
tors in the presence of SO interactions was initiated a 
rather long time ago^&S it continues to pose interesting 
and challenging problems. The Rashba SO interaction 6 
has received special attention, in part because of a pro- 
posal by Datta and DasiS that it could be exploited in 
a spin transistor — a device in which currents are mod- 
ulated by using a gate to alter the Rashba interaction 
strength. Some interesting refinements of the original 
idea have appeared in recent literature^ More gener- 
ally, there has recently been substantial theoretical work 
on spin-dependent transport in a 2DEG with the Rashba 
and^tl^r^yp^s^f^O intoactkin^ se e, e.g. Refs. p^ll^ . 



tions valid for weak SO interactions, which capture effects 
of the Rashba spin precession beyond the Dyakonov-Perel 
theory, have been derived and studiediiiiLi2iiS 

In this paper we derive a set of equations that pro- 
vide a complete description of coupled spin and charge 
diffusive transport in a 2DEG with the Rashba SO in- 
teractions. These equations capture a number of inter- 
related effects including spin accumulation and diffusion, 
Dyakonov-Perel spin relaxation and magnetoelectric and 
spin-galvanic effects.— This unified description is essen- 
tial, since spin transport is most easily detected in prac- 
tice through the spin accumulation it induces at the edges 
of the sample. 28 A complete understanding of the interre- 
lated spin accumulation and magnetoelectric effects in a 
given experimental situation can be obtained by solving 
the equations derived below. Our derivation is based on a 
microscopic evaluation of the disorder-averaged density- 
matrix response function, followed by an analysis of its 
long- wavelength, low frequency limit. We apply our 
equations to a simple model of electrical spin injection 
into a 2DEG from ferromagnetic spin-polarized contacts, 
placed on top of the 2DEG. We find that a voltage de- 
velops between ferromagnetic contacts when a spin po- 
larized current is injected into the 2DEG, that depends 
on the relative magnetization orientation of the contacts. 
Unlike all other known magnetoresistive effects in spin- 
tronics, this voltage drop is present even when the dis- 
tance between the electrodes exceeds the spin-diffusion 
length. 

The paper is organized as follows. In section[n]we out- 
line the density matrix response function formalism that 
we use to derive our spin and charge transport equations. 
In section lHll we derive transport equations from the low 
frequency, long-wavelength limit of the density matrix 
response function. Section llVl is devoted to a discussion 
of the relationship between formally exact equations of 
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motion for charge and spin-density operators in systems 
with Rashba spin-orbit interactions and arbitrary scalar 
disorder potentials and the coarse-grained dynamics pre- 
dicted by our diffusive transport equations. Section 
comments on the physical content of these equations and 
discusses an application to the case of electrical spin- 
injection from spin polarized contacts. Finally, in Section 
I VII we briefly summarize our findings. 



II. DENSITY-MATRIX RESPONSE FUNCTION 

Our analysis of coupled spin and charge transport in 
a semiconductor 2DEG system uses a model of noninter- 
acting electrons described by an effective-mass Hamilto- 
nian, moving in a random short-range spin-independent 
impurity potential. Because of the externally control- 
lable inversion-asymmetry of the quantum well confining 
potential, electrons experience a tunable SO interaction 
that we assume to be of the Rashba typepi The system 
is therefore described by a single-particle Hamiltonian 
H = Hq + Hi where 



# o = V ^ 



kcrcr' 



2m 



Xz 



x k] 



„t 



(1) 



is the effective-mass Hamiltonian with an additional 
Rashba SO interaction term. (We will use % = 1 units for 
convenience.) This interaction can be interpreted as Zee- 
man coupling to a k-dependent effective magnetic field 
2X(z x k). The impurity term in the Hamiltonian, 

Hi = j J2Vi(rM(r)**(r) = y £ ^(k-k')cL^. 

, ' r cr kk'er 

(2) 

describes the interaction of electrons with an impurity 
potential K(r) = u n J2 a <K r ~ r a)- The spin-independent 
random potential influences the electronic spin state by 
inducing transitions between momentum states that have 
different Rashba effective fields. The SO interaction lifts 
the spin degeneracy of the effective-mass Hamiltonian re- 
sulting in a momentum-dependent spin-splitting of the 
conduction band: 



e±(k) 



2m 



± Xk — /i. 



(3) 



We assume here that the Rashba spin-splitting is small 
compared to the Fermi energy Xkp <C £fj a good approx- 
imation in almost all cases of interest. 

Our analysis is based on an evaluation of the density- 
matrix response function using standard perturbation- 
theory methods^ The fundamental object in this ap- 
proach is the imaginary time Green's function 

Q aa . (r - r', r - r') = (T* a (r, r)**, (r', r')>, (4) 

where the angular brackets denote quantum, thermal and 
disorder averages. We compute the disorder averaged 



Green's function in the first Born approximation, which 
implies a self-energy in the Matsubara frequency repre- 
sentation given by 



= -7 



d 2 k 



(5) 



where 7 = n,Ug, n, is the density of impurities and Q± 
is the Green's function of the clean system without im- 
purities. The self-energy turns out to be momentum- 
independent and upon analytic continuation, iuj — > u) + 
irj, we obtain the familiar expression for the Born- 
approximation retarded self-energy: 



E CTCT / (ui + irj) = -—5 aa > 

ZT 



(0) 



where r = is the mean scattering time and go = 

m/n is the total density of states at the Fermi energy. 

It is convenient to decompose the disorder-averaged 
retarded and advanced real-time Green's functions G^'Jt 
into spin-independent singlet and spin-dependent triplet 
parts: 



Gf' A (k,u;)--r CTCT ,, (7) 



where r a are the usual spin-i Pauli matrices. The singlet 
and triplet Green's functions are given by 



Gf' A (k,c) = 



1 



w - £k - Afc ± £ u - & + Xk ± £ 



, (8) 



and 



G t (k,w) = — — 



u - 6c - Afc ± £ u - 6c + Afc ± ip 



(9) 



where 6c = k 2 /2m — /i. 

We can now proceed to evaluate the coupled spin and 
charge density response functions. We introduce the gen- 
eralized density operator 



£ CTlCT2 (r,t) = #t 2 (r,i)* CTi (r,i), 



(10) 



whose expectation value is the density matrix. (It is the 
matrix character of this quantity in spin-space that allows 
us to look at coupled spin-charge response; for present 
purposes it is adequate to specialize to diagonal elements 
in position-space.) From standard linear-response theory 
the retarded density response function is given by: 



X 



(T\(T2 -CT3 (7 A 
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,t-t') 



>',*')]>■ (11) 
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It is well known22, that this quantity can be evaluated 
to leading order in 1/kpt by summing all Born approx- 
imation self-energy and ladder vertex corrections to the 
polarization bubble (here I — vft is the mean- free-path). 
For (^-function impurities, the ladder sum for the Fourier- 
transformed retarded response function is a matrix geo- 
metric series which is easy to evaluate. We find that 

X(TlO'2,cr30'4( C lj ^) = 



simplifies to: 



1 



:f?0<5o-irT 3 8 a 



(12) 



where 



(q,fi) 

f d 2 k 



{2-Kf 



G£ ffl (k,0)G« 



and 



V =[1-1]- 



L (k + q,0), (13) 



(14) 



is the coupled spin-charge diffusion propagator or "diffu- 
son" . Summation over repeating spin indices is implied. 

At this point it is convenient to go to a more physical 
charge-spin-component representation for the response 
function: 



Xa/j(q,^) 



1 



(q.n)C. 



(15) 



where a, /3 = c, x, y, z. Inserting the identity matrix reso- 
lution ir? _ r Q , , between factors in Eg. 1121) . we obtain 

2 <*\&1 &2 a l 



Xapfaty = - — Y^ Icyfa ^)^7/3(q> fy~^QoS ai3 . (16) 

The integral over momentum in Ea. l|13fl is elementary 
but leads to cumbersome expressions for J(q, Q) that are 
listed in the appendix. We focus on the long- wavelength, 
low-frequency limit of the response function (|lfc>|) in the 
remaining sections of this paper. 



III. SPIN TRANSPORT EQUATIONS 

We are interested in the coupled dynamics of spin 
and charge, coarse-grained over lengths long compared 
to the mean-free path I and times long compared to the 
scattering time r. We concentrate here on the limit of 
weak SO interactions, A/cft <C 1 in which the scatter- 
ing time is much shorter than the spin-precession pe- 
riod and the spin-split Rashba bands are therefore not 
established. (The low-frequency, long-wavelength expan- 
sion of /(q, ft) is not analytic in the strong SO scattering 
limit.) In this diffusive limit, the inverse density fluctu- 
ation propagator (the diffuson) 2? _1 (q, £7) = 1 — /(q, 57) 



X> _1 (q,fi) = (-ifi 

/ iT sc q y 
ir sc q y 1/t_l 
-iT sc q x 

\ iT ss q x 




(17) 



where D = v f t/2 is the diffusion constant, t± = 
2t I (2Xkpr) 2 and t z = t±/2 are the in plane and out of 
plane spin relaxation times, and T sc — — 2A(Afci?r) 2 and 
r ss = AXeF t are the spin-charge and in-plane to out-of- 
plane spin couplings that result from the SO interactions. 
Note that iD/r^T 2 ^ = 1. 

Transforming this diffusion propagator to real space 
and time leads to the following system of coupled spin 
and charge transport equations that is the principal re- 
sult of this paper: 



— = DV 2 (N + g V c ) + 2T sc (z x V) • (S - g Q h) + I c , 



dN 

~di 

dS a 
dt 



= DV A - 



(S a - Q h a ) 



T ss p x V) x (S - g h)] a 

y(zxV) a (JV + ^c)+/ S '". 



(18) 



In these equations V c and h are the charge and spin (Zee- 
man) components of the external potential. The last 
term on the right-hand side of each equation has been in- 
serted by hand to represent charge and spin currents, I c 
and I s,a , vertically injected into the 2DEG. The factors 
of 2 and 1/2, that appear in front of the coefficient r sc , 
follow from the relationships between spin and charge 
densities and the corresponding combinations of elements 
of the density matrix. Note that in a generalization of 
the familiar Einstein relations, the external charge and 
spin potentials and the corresponding chemical poten- 
tials, N/ go and S/ go, are always summed; the charge and 
spin-densities respond as usual to electrochemical poten- 
tials and their gradients. (A 2DEG system with excess 
spin and charge densities N and S, has excess chemical 
potential (N ± 2|S|)/go f° r spins oriented along and in 
opposition to S respectively.) 

A physical understanding of the numerical values and 
the parametric dependences of the coefficients that ap- 
pear in front of the various terms in Eas. (|18|) is most 
easily obtained by considering the limit in which exter- 
nal potentials are absent. Then the drift and diffusion 
of charge and spin can be understood by considering 
the time evolution of electrons that start at the origin 
in specified spin-states and are scattered randomly be- 
tween various Rashba states at arbitrary angles on the 
Fermi circle. These electrons undertake random walks 
that make correlated steps of size ~ XkpT in spin-space 
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and I in position-space. The joint probability distribu- 
tion function that results from these correlated changes 
in spin and position is readily evaluated. Associating the 
coarse-grained spin and charge distributions with the dis- 
tribution of starting positions and spin orientations, the 
coefficients of n'th derivative terms in Eas. (|18(l arise from 
n'th order spatial moments of the spin and charge dif- 
fusion clouds. For example the diffusion constant D is 
related, as usual, to the second spatial moment of charge 
diffusion cloud and is therefore proportional to the square 
of the spatial step length I times the step rate r . Sim- 
ilarly T ss is due to spin-precession and is proportional 
to the first spatial moment of the S x spin projection in 
the diffusion cloud generated by spins that start with 
an orientation out of the plane. It is therefore propor- 
tional to the product of the spin-space and orbital-space 
step lengths and to the step rate. All non-standard co- 
efficients in our equations can be understood in terms 
of the correlation between velocity and spin-precession 
axis that exists throughout the random walk. This line 
of argument can be followed to provide an independent 
confirmation of Eas. (|18|l . 



IV. OPERATOR EQUATIONS OF MOTION 

Some insight into our general equations for the diffu- 
sive charge and spin density dynamics of 2DEG's with the 
Rashba spin-orbit coupling can be obtained by compar- 
ing Eas. (|I8() with the equations of motion of the charge 
and spin-density operators for this system. Let us first 
consider the Heisenberg equation of motion for the charge 
(or, more precisely, particle number) density operator 



N(t) = 

The equation of motion reads 
ON 



Of 



= i[H,N]. 



(19) 



(20) 



Since the particle number is conserved, we expect this 
equation to have the form of a continuity equation: 



ON 
~dt 



V • J c 



(21) 



where J c is the charge current density. Fourier trans- 
forming the charge density operator and evaluating the 
elementary commutator in Eq.J3U|) implies the following 
expression for the charge current density: 

J c = ~ (^V* CT - h.c.) + 2A (S x y - S*x) , (22) 



where 



S°(r) = -*t (r)r a CT ^ CT , (r)) 



(23) 



is the a-component of the spin density. The first term 
in Ea. (|22|l is the usual quantum- mechanical expression 



for the particle current density. We will call this con- 
tribution to the charge current a kinetic contribution. 
As seen from Eq.(J22J), SO interactions result in an addi- 
tional contribution to the charge current density, that we 
accordingly refer to as the spin-orbit contribution. This 
contribution is proportional to the in-plane spin densities. 
Comparing Eas. (|2T|l and J22J with the first of Eas. lfTHj) . 
we conclude that the kinetic contribution to the charge 
current transforms in the diffusive limit to a kinetic con- 
tribution of the standard form, proportional to both the 
diffusion constant and the electrochemical potential gra- 
dient. The SO contribution apparently remains separate 
and proportional to the spin-density, rather than being 
subsumed at long-wavelengths in the diffusive term. The 
total charge current is therefore given by, 



-DV(N 

2T SC [(S x 



QoVc) 
Qoh x )y 



(Sy - Qo h«)x] . (24) 



Let us now repeat the same analysis for the spin den- 
sity operators. Since the spin is not conserved, there is 
some freedom of choice in how the spin current density 
operator is defined^ We choose to define the spin cur- 
rent as the symmetrized product of the charge current 
discussed above and the spin operator, a definition that 
seems natural from a microscopic point of view and has 
been used previously, for example in discussing the spin- 
Hall effect M> The spin current density operator is, there- 
fore, also a sum of a kinetic and a spin-orbit contribution 
and has the following form: 



where 



kin 



3 S ' X = 


J kin < 




3 s ' y = 


K' y - 

kin 




3 S ' Z = 


-rS.Z 

**kini 




-—{ 

4m v 




,T a , - 
■' ' era' 



h.c. 



(25) 



(26) 



is the kinetic contribution to the spin current. Note that 
the current of the z-component of the spin has only a 
kinetic component. The Heisenberg equations of motion 
for the spin density operators can then be written in the 
following form: 



dS a 



-V-J s ' a + F a , (27) 
where F a is an additional source term that is given by: 



F. 



x.y 

F z 



= -2\mJ s x %, 
= 2\m(j s < x - 



Jy' V ) 



(28) 



As before, comparing Eqs. (J2SI>— (EHJ with Eq.(UH|), we 

conclude that the kinetic contribution to the spin cur- 
rents is proportional to the gradient of the spin electro- 
chemical potential and, in addition, that the currents of 
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the in-plane spin components have SO contributions as 
in the microscopic equations of motion: 



JS ,x = _ DV ( S *_ eoh *) + _!Z.( N + eQVe )j jt 



J'rv = -DV (S v - Qohv) - -^(JV + Qo V c )x, 
J s ' z = -DV(S Z - g h z ). 



(29) 



As in the charge current case, the spin-orbit contribu- 
tion to the microscopic spin current is not subsumed in 
the diffusive contribution, but appears separately. In- 
terestingly, a change occurs in passing from the micro- 
scopic expression to the coarse-grained transport theory 
expression in that the chemical potential (proportional 
to the density N in 2D) is replaced by the electrochemial 
potential. Because the spin-orbit spin current contribu- 
tions are proportional to the charge density in the ab- 
sence of external fields, they are non-zero in equilibrium, 
as noted by Rashba. 36 Although the constant equilibrium 
spin currents in a uniform system have no physical con- 
sequences as far as we are aware, these spin-orbit terms 
in the spin current do play an important role in cou- 
pled spin-charge transport as we illustrate in the follow- 
ing section. Also note that the source terms F a in the 
microscopic equations of motion appear in an almost, but 
not completely, identical way in the diffusive equations of 
motions; the constants multiplying the currents are twice 
as large in the diffusive case: T ss /D = 4Am instead of 
2Am in Ea. (|28|l and N is replaced by its electrochemical 
equivalent N + QqV c . 

Finally, let us comment on the relationship be- 
tween our results and the recently discovered spin-Hall 
effect, 2fi i 31 &&mMm3MML i.e. a transverse J s ' z cur- 
rent in response to an in-plane electric field. Note that, 
as in the microscopic expressions, the current of the z- 
component of the spin has only a kinetic contribution. 
This means that, apparently the spin-Hall effect does 
not occur in the diffusive limit, since the current of S z 
in Ea. pHfl has only a diffusive contribution, that does 
not react to the electric field. In the diffusive regime, 
the spin current divergence (from the spin-orbit contri- 
bution), produced by a uniform electric field, is balanced 
in the steady-state by spin relaxation. This balancing 
leads to a spin-polarization perpendicular to the electric 
field direction, as noted some time agoii 2 *^ In the regime 
of resolved spin-orbit induced spin-splitting, the diffusive 
transport picture of Ea. i|18|) is no longer applicable. In 
this case, the current of the ^-component of the spin will 
have a contribution, proportional to z x E, due to the 
intrinsic spin-Hall effect. 



V. APPLICATIONS OF COUPLED 
SPIN-CHARGE TRANSPORT EQUATIONS 

In the case of an infinite 2DEG , Eqs.JTSJ can be solved 
by Fourier transformation. Rotating coordinate axes so 



that the y-axis is along the direction of q, brings the in- 
verse diffusion propagator Ea. (|I7|) to the block-diagonal 
form: 

£> _1 (q, fl) = (-ifl + Dq 2 ) 1 + 

iT sc q 

iY sc q 1/tx 

1/t x -iT ss q 

iT ss q 1/t z 



(30) 



The eigenmodes are then easily calculated to be: 




The ifli± modes correspond to coupled diffusion of 
charge and the in-plane spin density component that is 
transverse to the direction of q, i.e. S x in this convention. 
Note, that the mode if2i_ is gapless at q = 0. This means 
that this mode corresponds to a (nearly) conserved quan- 
tity, with a very long relaxation time at small q. Exactly 
at q = this quantity is of course simply the conserved 
total particle number. However, at finite wavevectors it 
corresponds to a linear combination of the charge density 
and the x-component of the spin density. 

The iU,2± modes correspond to coupled diffusion of 
S y and S z spin densities. This coupling originates from 
the Rashba spin precession as explained above. Note 
that iU,2- has a minimum at a finite wavevector q* — 
\/l5Am/2, as discovered previously in Ref.0]. This 
means that the S v ' z Fourier component with the slow- 
est relaxation rate will actually be at q = q* , unlike in 
the case of the ordinary diffusive relaxation, where the 
slowest relaxation rate is at q = 0. 

Let us now look at stationary solutions of Eq. Ijf8|l . For 
simplicity and clarity of presentation we will assume that 
spin and charge densities are uniform in the ^-direction, 
and, therefore, the inverse diffusion propagator has the 
simple block-diagonal form Eq. l|3l)[) . In the following 
paragraphs we discuss the stationary state response of 
the 2DEG system to external spin and charge currents 
injected or drained along lines of constant y as illustrated 
in Fig. We first consider the response to a flux of the 
z-component of the spin, l s,z . Inverting the lower block 
of the inverse diffusion propagator Ea. (|3U[l . we obtain: 



S z (q) = I s 
and 
S y (q) = r 



Dq 2 + 1/tx 



D 2 q 4 - 4(Ae F r) 2 q 2 + 32(X 2 me F T) 2 ' 

■ -iTssq 

L> V - 4(Ae F r) 2 q 2 + 32(A 2 me F r) 2 ' 



(32) 



(33) 
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injector 



detector 



transport equations: 




FIG. 1: Cartoon of the spin injection experiment. A po- 
larized current, that has a polarization component along the 
x-direction (i.e. along the electrode) is injected at the left 
electrode and collected at the right one. A voltage develops 
between the electrodes that depends on the spin-polarizations 
of both emitter and collector. 



The Fourier transform to real space is readily evaluated 
by contour integration. Poles occur at the roots of the 
denominator located at 



q = ±Xm J 2(1 ±i\fl). 



(34) 



Note that all the roots are complex. This means that the 
nonequilibrium spin density profile in this case will not 
have the usual form, exponentially decaying away from 
the point where the current is injected, with a character- 
istic spin diffusion decay length. (Because the spin-orbit 
coupling is linear in momentum, the distance travelled 
by a Fermi energy electron during one spin-precession, 
oc 1/Am, is independent of the Fermi momentum.) In- 
stead it will clearly involve an oscillatory component, 
which is the remnant of the Rashba spin precession in 
the diffusive regime. The inverse of the characteristic 
decay length and the inverse period of the spatial os- 
cillations are given by the imaginary and real parts of 

\m\J 2(1 + correspondingly. 

Let us now turn to the more interesting issue of signa- 
tures of the spin-charge coupling in our transport equa- 
tions in spin injection experiments. We imagine the ge- 
ometry schematically depicted in Fig. ^ Assume two 
infinitely long ferromagnetic electrodes are placed on top 
of the 2DEG sample a distance L from each other. Let a 
charge current /, polarized in the x-direction (i.e. along 
the electrode), be injected into the 2DEG from the y = 
electrode. Assume that the degree of spin polarization 
of this current is a, i.e. the injected spin current is 
js,x _ a j Assume in addition that this current is ex- 
tracted at the second electrode at y = L, which has a 
degree and sign of spin polarization denoted by /3, that 
can differ from a. This circumstance is expressed com- 
pactly by the following source terms in our spin-charge 



r 

I S,X 



(y) = I[S(y)-S(y-L)}, 
(y) = I[a5(y)-06(y-L)]. 



(35) 



We now evaluate voltage, i.e. the electrochemical poten- 
tial difference, that develops between the two electrodes 
and discuss how it depends on the relative spin polariza- 
tion of the electrodes. 

Inverting the upper block of the diffusion propagator- 
Ed. (I30|) . we obtain the Fourier transformed local elec- 
trochemical potential change U = N/ go + V c that is 
generated in the 2DEG in response to the injected spin- 
polarized current: 



QoU(q) 



II 



-iqL\ 



1 



Dq 2 



- l(a~f3e- zqL ) 



2tT, 



D 2 q {q 2 + 1/Dt±) ' 



(36) 



We define the effective electric field in the 2DEG in the 
usual way in terms of the gradient of the electrochemical 
potential: 



E= --VU. 

e 



(37) 



The electric field response to the injected current can 
then be easily calculated by an inverse Fourier transfor- 
mation: 



E (y) = -^j d qq u( q y 



I 



eg D 

(38) 

The voltage between the electrodes is therefore given by: 

r-L 



V = 



IL 

eg D 



" 



E(y)dy 



(39) 



Note that in addition to the usual Ohm's law contribu- 
tion to the potential drop, the first term on the right 
hand side of Eq. l|39|l . there is also a contribution propor- 
tional to the difference of spin polarizations of the two 
ferromagnetic electrodes. This sensitivity of the resis- 
tance of a paramagnetic system to the spin-polarization 
of the current-carrying electrons results from the spin- 
charge coupling terms in our transport equations that 
generate a contribution to the electric field proportional 
to r sc (see Eg. (13811 ). Note that this contribution to the 
voltage is present even when the electrodes are separated 
by a distance larger than the spin-coherence length. 

This effect could be studied by attaching voltage 
probes to the 2DEG near the ferromagnetic electrodes, 
or simply by measuring the voltage drop between the 
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ferromagnetic electrodes. In the latter case, the total 
voltage will contain contact contributions from the local 
voltage drops between the ferromagnetic electrodes and 
the 2DEG. The spin-orientation dependent voltages that 
we discuss will, in general, need to be distinguished from 
other spin polarization-dependent voltages that occur in 
magnetotransport, for example the spin-polarization de- 
pendent open-circuit voltages measured by ferromagnetic 
electrodes first discovered in the seminal work of John- 
son and SilsbeeiSI The effects that we discuss here can 
be distinguished in several ways. First of all, the voltage 
differences that we have calculated are ones that would 
be measured by paramagnetic voltage probes. Secondly, 
our voltages have a characteristic dependence on the po- 
larization of the injected spin current. In the case of the 
effect described by Ea. (|39[l . the voltage drop is maximal 
when the first electrode is polarized along the x-direction, 
while the second one is polarized along the — x-direction. 
On the other hand, if the electrodes are polarized along 
the y and — y-directions, the polarization-dependent volt- 
age drop will vanish. 



VI. CONCLUSIONS 

In this paper we have examined the issue of the diffu- 
sive spin and charge density transport in Rashba 2DEG 
systems. The separation between the spin and momen- 
tum relaxation time scales in the diffusive regime has 
allowed us to use a statistical description, where the spin 
and charge transport is described by local spin and charge 
electrochemical potentials and their gradients. Our the- 
ory thus generalizes the usual two-component theory of 
diffusive spin transport that has found numerous 
successfull applications, in particular in the theory of 
spin-dependent transport in magnetic multilayers. 43 

Our equations with appropriate boundary conditions 
can be used to model experiments on coupled spin-charge 
transport in 2DEG systems with the Rashba SO interac- 
tions, involving both electrical and optical spin injection. 
By comparing our equations, valid in the diffusive trans- 
port regime, with the exact operator equations of motion, 
we have inferred relationships between spin and charge 
current densities and spin and charge electrochemical po- 
tentials and their gradients. These expressions can be 
used to devise appropriate boundary conditions that are 
necessary to supplement our transport equations at the 
2DEG boundaries. 

As an example of the application of our equations, we 
have considered a simple electrical spin injection experi- 
ment and shown that a voltage will develop between two 
ferromagnetic contacts if a spin-polarized current is in- 
jected into a 2DEG sample, that depends on the relative 
magnetization orientation of the two contacts. Unlike the 
giant magnetoresistance and other familiar magnetore- 
sistive effects in spintronics, this voltage drop is present 
even when the distance between the electrodes exceeds 
the spin-diffusion length. 
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VII. APPENDIX 

This appendix summarizes some technical details of 
the density matrix response function calculation in sec- 
tion [H] The main technical problem is the evaluation of 
the matrix elements /, 



0102,0304 



(q, which turn out to 



have the following general form: 

7xi02, 03,04 (Qi ^) — I (Q; ^)^(TiCT3 ^CT2CT4 

/«(q,0)+/«(q,fi) 



-T q T q 
0-30-1 'o- 2 o 4 



/**(q,Q)-7«(q,fi) „ 



where 



and 



-up 



-ie- % f 



Hi 







(40) 



(41) 



(42) 



are the Pauli matrix components along and perpendicu- 
lar to the direction of q (ip is the angle between q and 
the x-axis). The explicit expressions for the functions 
7 SS ,/ S *,7" and J", that appear in the matrix elements 
of J, are: 



r s ( q ,n) 



2 1 1 

/o /- /+ 



(43) 



7 st (q,0) 



A v / 2Dq 2 r 



/ i 2L>q 2 r / 2Dq 2 T 



iq\T{\kpT) 2 , 
(44) 



J?(cbfi) = 4 



2 1 



1 



(45) 
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where 



J«(q,fi) 
2 



1 



4L>q 2 r 

f 2 - Dq 2 r f_ - Dq 2 r f\ - Dq 2 r 



f- 



fo 

f± 



y/(l-iriT) 2 + 2Dq 2 T, 

V(l - ifir ± ifl so T) 2 + 2Dq 2 r, 



(46) 



(47) 



and Q so = 2\kp is the Larmor precession frequency asso- 
ciated with the Rashba field. The spin-charge coupling is 



generated by the last term in the expression for I st (q, 57). 
Note that, unlike other elements of the matrix /(q, 57), 
the term responsible for the spin-charge coupling can only 
be calculated perturbatively in A and q. 



Expanding Eas. l|43(l - (|47|) to leading order in itlr, 
Dq 2 r and XUft, one arrives at Ea. (|18|l . which we have 
concentrated on in this paper. These expressions should 
be useful in their original form, however, in systems with 
the Rashba interactions strong enough that the spin- 
splitting of the band energies is not smeared out by dis- 
order. 
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